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ABSTRACT 

We report on the observations of 14 dwarf spheroidal galaxies with the Fermi Gamma-Ray 
Space Telescope taken during the first 1 1 months of survey mode operations. The Fermi tele- 
scope, which is conducting an all-sky y-ray survey in the 20 MeV to >300 GeV energy range, 
provides a new opportunity to test particle dark matter models through the expected y-ray emis- 
sion produced by pair annihilation of weakly interacting massive particles (WIMPs). Local 
Group dwarf spheroidal galaxies, the largest galactic substructures predicted by the cold dark 
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matter scenario, are attractive targets for such indirect searches for dark matter because they are 
nearby and among the most extreme dark matter dominated environments. No significant y-ray 
emission was detected above 100 MeV from the candidate dwarf galaxies. We determine upper 
limits to the y-ray flux assuming both power-law spectra and representative spectra from WIMP 
annihilation. The resulting integral flux above 100 MeV is constrained to be at a level below 
around 10~^ photons cm~^s~^Using recent stellar kinematic data, the y-ray flux limits are com- 
bined with improved determinations of the dark matter density profile in 8 of the 14 candidate 
dwarfs to place limits on the pair annihilation cross-section of WIMPs in several widely studied 
extensions of the standard model, including its supersymmetric extension and other models that 
received recent attention. With the present data, we are able to rule out large parts of the pa- 
rameter space where the thermal relic density is below the observed cosmological dark matter 
density and WIMPs (neutralinos here) are dominantly produced non-thermally, e.g. in models 
where supersymmetry breaking occurs via anomaly mediation. The y-ray limits presented here 
also constrain some WIMP models proposed to explain the Fermi and PAMELA e^e~ data, 
including low-mass wino-like neutralinos and models with TeV masses pair-annihilating into 
muon-antimuon pairs. 



1. Introduction 



A wealth of experimental evidence and theoretical arguments have accumulated in recent years in favor 
of the existence of some form of non-baryonic cold dark matter (CDM) to explain the observed large-scale 
structure in the Universe. According to the most recent estimates, CDM constitutes approximately one- 
fourth of the total energy density of the Universe. However, very little is known about the underlying nature 
of this dark matter, despite the efforts of high-energy physicists, astrophysicists and cosmologists over many 
years, and it remains one of the most fascinating and intriguing issues in present day physics. One appealing 
possibility is that CDM consists of a new type of weakly interacting massive particle (WIMPs), that are 
predicted to exist in several theories beyond the Standard Model of particle physics. Such WIMPs typically 
have pair-annihilation cross sections that, for their natural mass range (between a few GeV and a few TeV), 
drive a thermal relic abundance in the same ballpark as the inferred amount of cosmological dark matter. 
Pair-annihilation (or decay) would also occur today, yielding, among other particle debris like energetic 
neutrinos, (anti-)protons and electron-positron pairs, a significant flux of high-energy y-rays. If the dark 
matter is meta-stable, its decay products would also produce potentially detectable y rays. While the results 
we present here would constrain this type of scenario, we assume here that the dark matter particle is stable. 

Cosmological N-body simulations of structure formation show that the dark matter halos formed by 
WIMPs are not smooth and have a large number of bound substructures (subhalos) who se numbers increase 
with decreasing mass dSpringel et al.ll2005l : iKuhlen et al.ll2008l : iDiemand et al.ll2005h . Dwarf spheroidal 
galaxies (dSphs), the largest galactic substructures predicted by the CDM scenario, are ideal laboratories for 
indirect searches for dark matter, through the observation of dark matter annihilation (or decay) products, 
for the following reasons. The mass-to-light ratios in dSphs can be of order 100 - 1000 (see Table [D, 
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showing that they are largely dark matter dominated systems. Therefore, the stars serve as tracer particles in 
the dark matter gravitational potential and the dark matter distribution in these dwarfs may be constrained 
using stellar kinematics. In addition, dSphs are expected to be relatively free from y-ray emission from other 
astrophy sical sources as they have no detec t ed neutral or ionized gas, and little or no recent star formation 
activity (IMatedllQQSuGallagher et al.ll2003l : iGrcevich & Putmanl l2009). thus simplifying the interpretation 
of a gamma-ray e xcess that would b e detected in the direction of a dSph. In addition, the Sloan Digital Sky 
Survey (SDSS, see I York et al has led, in recent years, to the discovery of a new populati on of ultrafaint 



Milky Way satellites, comprising about as many (new) objects as were previously known (|Willman et al. 



2005 



Zucker et al.ll2006l : llrwin et al.ll2007l : IWalsh et al.ll2007l : lBelokurov et al.ll2007h . This new population 



of extremely low-luminosity, but dark matter dom inated galaxies could in particular be very interesting for 
indirect dark matter searches (IStrigari et al. I I2OO8I) . especially with the Fermi Gamma-Ray Space Telescope. 



Fermi is a new generation space observatory, which was successfully launched on June 1 1th, 2008, and 
has been operating in nominal configuration for scientific data taking since early August 200qj. Its main 
instrument, the Large Area Telescope (LAT), is designed to explore the high-energy y-ray sky in the 20 MeV 
to > 300 GeV energy range, with unprecedented angular resolution and sensitivity. Several studies have been 
performed to determine the sensitivity of the LAT to dark matter annihilation signals (e.g. iBaltz et al.ll2008h 
and the FermZ-LAT collaboration is currently exploring several potentially complementary search strategy 
for y-ray emission from dark matter. We focus here on the search for a y-ray signal in the direction of a 
selection of 14 dSphs. The criteria for this selection, together with the description of the Ferm/-LAT data 
analysis, are presented in ^ We determine upper limits to the y-ray flux employing both power-law spectra 
with spectral indexes in the range between 1 and 2.4 ( 92. lb . and spectra resulting from the annihilation of 
several representative WIMP models, for various WIMP masses ( 93.21 ). To turn these results into limits on 
the WIMP pair annihilation cross section, we restrict our focus to a subset of 8 dSphs that are associated 
with stellar data of good enough quality to allow for an accurate modeling of their dark matter content. 
We then present an updated determination of the assumed Navarro-Frenk- White dark matter density profile, 
using a Bayesian analysis and up-to-date stellar velocity data ( 93.11 1. We show in 93.21 results for WIMP 
models in the context of minimal supergravity, of a general weak-scale parameterization of the minimal 
supersymmetric standard model, of the minimal anomaly-mediated supersymmetry breaking scenario, and 
of universal extra-dimensions. Finally, in 93.31 we discuss these constraints in the context of dark matter 
annihilation models that fit the PAMELA and Fermi e^e" data, putting special emphasis on the effect that 
a possible contribution of secondary y-ray emission from IC scattering has on them. Our main conclusions 
are summarized in 94] 



2. Fermi-LXT observations and data analysis 

The LAT i s an electron-positro n pair conversion telescope sensitive to photon energies from 20 MeV 



to > 300 GeV (lAtwood et al. 



20091) . It is made of 16 towers each comprising a tracker and a calorimeter 



' For more details, see the Fermi website at |http://ferrm.gsfc.nasa.gov71 
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Table 1 . Properties of the dwarf spheroidals used in this study. 



Name Distance year of M1/2/L1/2 1 b Ref. 

(kpc) discovery ref. 8 



Ursa Major II 


30± 5 


2006 


4000+""" 

4UUU 2,00 


152.46 


37.44 


1,2 


Segue 2 


35 


2009 


650 


149.4 


-38.01 


3 


Willman 1 


38± 7 


2004 


' '"-440 


158.57 


56.78 


1 


Coma Berenices 


44± 4 


2006 


11 00+™ 


241.9 


83.6 


1,2 


Bootes II 


46 


2007 


18000?? 


353.69 


68.87 


6,7 


Bootes I 


62±3 


2006 


1 'OU_7QQ 


358.08 


69.62 


6 


Ursa Minor 


66± 3 


1954 




104.95 


44.80 


4,5 


Sculptor 


79± 4 


1937 


18^6 


287.15 


-83.16 


4,5 


Draco 


76± 5 


1954 




86.37 


34.72 


4,5,9 


Sextans 


86± 4 


1990 


120+^^ 


243.4 


42.2 


4,5 


Ursa Major I 


97±4 


2005 


1800!^300 


159.43 


54.41 


6 


Hercules 


132±12 


2006 


MOO^ig"" 


28.73 


36.87 


6 


Fornax 


138±8 


1938 


7+2.8 


237.1 


-65.7 


4,5 


Leo IV 


160±15 


2006 


260+'™ 


265.44 


56.51 


6 



Note. — M1/2/L1/2 is the ratio of the total mass within t he 3D half-Ught radius to the 
stellar luminosity within the same radius from IWolf et al] j2009l) . The problematic re- 
sult for Bootes II is further discussed in the text. Uncertainties in the determination of 
this mass-to-light ratio (unavailable for Bootes II and Segue 2) arise from the errors in 
both Ml/2 and L|/2, but they do not change the qualitative conclusion tha t these dSphs 
are da rk matter dominated even with in their stellar extent. R efe rences: (l)IStrigari et alj 



1200a 



<2008l). (2) Simon & Geh3 l2007h. OllBelokurov et al.l |2009|), (4 ) |Penarrubia et alj . . 
(S'tlMated il99A (6 ) lMartin et alj lioOSi) , (7) lKoch et alj j2009l) (8) IWolf et alj |2009|) (9) 
iBonanos etalJto04l) 
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undemeath. The tracker is made of silicon-strip active planes interleaved with tungsten foils, and is respon- 
sible for the conversion of the incident photon into an electron-positron pair and for the tracking of the latter 
charged particles. The energy of the photon is mainly estimated from the light deposited in the CsI(Tl) 
scintillators that constitute the calorimeter. Finally, an anticoincidence detector, made of more than 100 
plastic scintillators, covers the 16 towers in order to be able to reject the charged particle background. The 
LAT nominally operates in a scanning mode observing the whole sky every 3 hours, the resulting overall 
coverage of the sky being fairly uniform. The analysis described here uses data taken in this mode during 
the first eleven months of sky survey operation (August 4 2008 to July 4 2009). 

The dSphs that have been considered for this work are Usted in Table [T] They were selected based 
on their proximity, high Galactic latitude and their dark matter content, which have been estimated from 
the most recent resolved stellar velocity measurements. In particular. Carina and Sagittarius were discarded 
based on the fact that they are at low latitude {\b\ < 30°), and thus subject to potentially large systematics 
due to uncertainties associated with the modelin g of the Galactic D iffuse emission as seen by the LAT. In 



addition. Segue 1 is a controversial case : while iGeha et all (|2009ah concluded tha t this dwarf is the most 



promising satellite for indirect dark matter detection, this claim was challenged by iNiederste-Qstholt et al. 



who contend that Segue 1 is a star cluster stri pped early o n from the Sagittarius galaxy. As new 



stellar data on Segue 1 are currently being analyzed (iGehal 120091) . which may greatly improve the still 
uncertain measurements of the density profile, we defer its study to an upcoming dedicated publication. 
Finally, Bootes II is modeled based on the published data on 5 stars. The result that we obtain for the mass 
to luminosity ratio is unrealistic, and could mean that something is wrong with the stellar membership of 
this system or that it is simply unbound. Nevertheless, more data exist and are currently being reduced, and 
proposals are under way to increase the stellar dataset for this object. As a result, we keep Bootes II in the 
present analyis, since in the future it may prove to be one of the best candidate dwarfs. Finally, in section |3l 
we use a subset of 8 dwarfs that have robust stellar kinematic data to further constrain models. 

The dat a reduction ma k es use of the standard LAT ground processing and background rejection scheme 



described in lAtwood et all (120091) . and we consider only 'Diffuse' class events, which have the highest 
probability of being photons. Throughout, we use the Fermi ScienceTools version v9rl5, a software 
package dedicated to the Fermi-hAT data analysi^l. First, observations toward each dSph are extracted in 
14 regions of interest (ROI), by keeping events that have a reconstructed direction of incidence at most 10° 
away from the dwarf position. This ROI accomodates the large point spread function (PSF) of the LAT 
at low energy. Indeed, the LAT PSF, which depends on the photon energy and angle of incidence, can 
be approximated by the function 0.8 (S/lGeV)""-^ deg, yielding ~ 5" at 100 MeV. Next, in order to avoid 
calibration uncertainties at low energy and background contamination at high energy, we apply a cut on the 
reconstructed energy E: E > 100 MeV and E < 50 GeV. Here we employ a somewhat conservative cut 
at high energies to reduce the background, but work is ongoing within the Fermi collaboration to develop 
an improved event selection which will have less high energy background contamination (Abdo et al., in 
preparation). To avoid albedo y-ray contamination, we also select Good Time Intervals (GTIs) when the 



^http://fermi. gsfc.nasa.gov/ssc/data/analysis/software/ 
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entire ROI is above the albedo horizon of the Earth (105° below the zenith). Furthermore, the Earth limb 
appears at a zenith angle of 1 13° from Fermi's orbit. Thus, time periods during which the spacecraft rocking 
angle (the angle between LAT normal and Earth-spacecraft vector) is larger than 43" are excluded as an 
additional guard against Earth albedo y-ray contamination. 

The resulting dataset is analyzed with a binned likelihood technique fcashlll979l : iMattox et al.lll996h . 
implemented in the ScienceTools as the gtlike task, gtlike uses the maximum likelihood statistic to 
fit the data to a spatial and spectral source model. Because the numbers of photons from sources near the 
detection limits are fairly small, gtlike calculates a likelihood function based on the Poisson probability 
using the source model folded through the LAT Instrument Response Functions (IRFs) to provide the 
expected model counts. This analysis relies on version P6_V3 of the IRFs. For each ROI, the source 
model includes all the preliminary 11 month LAT catalogue sources within a 10° radius of each dwarf. 
Following the analysis procedure used in the development of the LAT catalogue, these sources are modeled 
as point-like with power-law spectra, a reasonable approximation in the absence of dedicated studies for 
most of them. Furthermore, the positions and spectral parameters of these sources are being kept fixed 
during the fitting procedure. It also includes the models currently advocated by the LAT collaboratior0 
for the Galactic diffuse emission and for the corresponding isotropic component (which accounts for any 
extragalactic diffuse emission and any residual charged background contamination). Their independent 
normalizations are kept free during the fit in order to account for uncertainties in modeling these diffuse 
components. Finally, the dSphs are modeled as point sources, their localization being kept fixed during the 
fit. Given the limited angular resolution of the LAT and the limited statistics, the point source approximation 
is reasonable for all the selected dwarfs. To model the source spectra, we employ two strategies: we employ 
model-independent power-law spectra, with a wide range of spectral indexes from 1 to 2.4, discussed in 
92.11 as well as a collection of motivated and representative y-ray spectra from WIMP pair annihilation, for 
wide ranges of masses and for several WIMP models (see 93.21 ). 



2.1. Power-law modeling 

For a power-law spectrum, the differential flux is written as 

dN ( E 

= A^J _ , (1) 

dEdAdt \Eol 

where E is the photon reconstructed energy in the restricted energy range 100 MeV to 50 GeV, and Eq is 
an arbitrary energy scale set to 100 MeV. Such a model has two unknown parameters, the photon index F 
and the normalization parameter Nq. In this analysis, we fix F to five possible values, F = 1, 1.8, 2.0, 2.2 
and 2.4. While the last four indices provide constraints on standard astrophysical source spectra, the very 

^http://www-glast.slac. stanford.edu/software/IS/glas t_lat_performance.htm 

''Detailed description can be found under 'Model Description' at the following web page : 
|http://fernii.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.htrnl| 
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hard index of F = 1 is motivated by the dark matter annihilation models in lEssig et alj (I2009I) . A^o is fitted 
to the data in each ROI separately, together with the isotropic and Galactic diffuse normalizations. The best 
fit values and corresponding errors of these three parameters, for the case T = 2, are gathered in Table |2l In 
the third column, the errors on A^o are several orders of magnitude larger than the fitted values, which means 
that the latter are compatible with zero. We thus conclude that no significant signals from the direction of 
the selected dSphs are detected. The same conclusion is reached for the other values of the photon index 
F. 

We expect variations across the ROIs around the ideal value of 1 for the normalizations of the two 
diffuse components, due to possibly slightly different exposure, statistical fluctuations, and inaccuracy of 
the underlying spatial model and spectral shape. Nevertheless, the fit results for the isotropic component 
remain close to 1, while the deviations for the Galactic diffuse component are somewhat larger, which is 
expected as spatial inaccuracies of the model are expected to vary from one ROI to another. Fig. [T] shows 
examples of the spectral fits to the data and the residuals of these fits for WilLman 1, which has the largest 
fit residuals, and Draco, which has residuals typical of most of the fits. 



Flux upper limits are then derived, based on a profile likelihood technique (lBartlettlll953l : iRolke et al. 



20051) . In this method, the normalization of the power law source representing the dwarf galaxy is scanned 
away from the fitted minimum while the remaining two free parameters (the normalizations of the two 
diffuse backgrounds) are fit at each step. The scanning proceeds until the difference of the logarithm of the 
likelihood function reaches 1.35, which corresponds to a one-sided 95% confidence level (C.L.). Extensive 
work to test this method, using Monte Carlo simulations as well as bootstrapping the real data, indicate 
that profile likelihood method as implemented in the Science Tools is slightly overcovering the expected 
confidence leveH As a result, we believe the ULs presented in this paper are conservative. 

In Table |3l we report flux upper limits in two different energy ranges (from lOOMeV to 50GeV and 
from IGeV to 50GeV), for the five different values of F. As expected, the Fen«/-LAT limits are stronger 
for a hard spectrum which predicts relatively more photons at the dwarf location at high energy. At higher 
energies, the LAT PSF is significantly reduced while the effective area is significantly higher; in addition, 
the diffuse backgrounds have relatively soft spectra compared to all but the softest models considered. For 
example, for a power law index of 1 when analyzing the full energy range, the upper limits are roughly 10 
times lower than for an index of 2.2. 

As can be seen from Table 3, the different dwarfs give roughly similar gamma-ray flux upper limits, as 
expected given the fairly uniform Fermi exposure across the sky. However, the limits do vary from dwarf to 
dwarf due to, for example, the direction dependence of the diffuse background intensity and the proximity 
of bright gamma-ray sources. In general, Ursa Minor gives the lowest flux limits while Sculptor, which is 
within a couple of degrees of a bright background point source, gives the highest flux limits. 



^The upper limits derived from the tests covered 98% of the trials instead of the 95% required. 



Table 2. Results for the fit of the three normalization parameters in each ROI. 



dSph 


Galactic foreground 


Isotropic component 










xur^) 






Ursa Major II 


0.89 ± 0.07 


0.97 ± 0.06 


5 


8 


X 


lo- 


09 


± 9.29 X 


10- 


-06 


Segue 2 


1.01 ±0.04 


1.03 ±0.06 


2 


3 


X 


ur 


15 


± 5.22 X 


10- 


-10 


Willman 1 


0.44 ± 0.25 


1.07 ±0.07 


7 


3 


X 


ur 


05 


± 1.21 X 


10- 


-04 


Coma Berenice 


0.90 ±0.15 


1.06 ±0.06 


2 





X 


10- 


13 


± 4.54 X 


10- 


-09 


Bootes II 


0.96 ±0.14 


1.19 ±0.08 


3 


4 


X 


10- 


12 


± 2.69 X 


10- 


-08 


Bootes I 


0.80 ±0.13 


1.26 ±0.08 


1 


2 


X 


10- 


12 


± 2.32 X 


10- 


-08 


Ursa Minor 


0.50 ±0.11 


1.13 ±0.06 


4 


2 


X 


10- 


13 


±6.15x 


10- 


-09 


Sculptor 


0.53 ±0.15 


1.03 ±0.06 


1 


2 


X 


10- 


04 


± 1.40 X 


10- 


-04 


Draco 


0.70 ± 0.09 


1.08 ± 0.06 


2 


1 


X 


10- 


11 


± 5.68 X 


10- 


-08 


Sextans 


1.00 ±0.10 


1.09 ±0.06 


6 


6 


X 


10- 


12 


± 3.15 X 


10- 


-08 


Ursa Major I 


0.71 ±0.27 


1.02 ±0.07 


2 


9 


X 


10- 


14 


± 1.84 X 


10- 


-09 


Hercules 


0.93 ± 0.07 


1.27 ±0.09 


1 


1 


X 


10- 


14 


± 2.56 X 


10- 


-09 


Fornax 


1.01 ±0.16 


0.86 ± 0.06 


1 


1 


X 


10- 


15 


±2.31 X 


10- 


-10 


Leo IV 


0.94 ±0.11 


1.23 ±0.06 


2 


1 


X 


10- 


11 


± 8.46 X 


10- 


-08 



Note. — Results are shown for the case T = 2. EiTors are statistical only. Nq is the power-law 
prefactor in eq.(T). 



Table 3. Flux Upper Limits at 95% C.L. 



£>100MeV £>lGeV 



Spectral index T 


1.0 


1 


.8 


2.0 


2.2 


2.4 


1.0 


1 


.8 


2.0 


2.2 


2.4 


Ursa Major II 





14 


1 


28 


2 


15 


3 


41 


5 


12 





09 





23 





29 





33 





37 


Segue 2 





10 





71 


1 


28 


2 


33 


4 


21 





06 





12 





15 





17 





19 


Willman 1 





14 


1 


63 


3 


02 


5 


22 


8 


39 





11 





34 





40 





44 





47 


Coma Berenices 





08 





43 





69 


1 


11 


1 


74 





05 





07 





08 





09 





09 


Bootes II 





13 





77 


1 


19 


1 


77 


2 


53 





08 





13 





14 





14 





15 


Bootes I 





12 


1 


02 


1 


71 


2 


70 


3 


96 





09 





23 





28 





31 





33 


Ursa Minor 





08 





39 





60 





88 


1 


26 





05 





08 





09 





10 





11 


Sculptor 





22 


2 


43 


3 


88 


5 


71 


7 


76 





12 





34 





39 





44 





48 


Draco 





09 





59 





94 


1 


41 


1 


94 





06 





13 





16 





21 





26 


Sextans 





09 





56 





97 


1 


62 


2 


55 





06 





10 





13 





16 





21 


Ursa Major I 





09 





48 





77 


1 


23 


1 


90 





06 





09 





10 





12 





14 


Hercules 





33 


1 


51 


2 


22 


3 


23 


4 


63 





24 





30 





30 





28 





27 


Fornax 





12 





94 


1 


72 


3 


05 


5 


04 





09 





14 





16 





17 





18 


Leo IV 





12 





96 


1 


58 


2 


47 


3 


64 





08 





21 





26 





32 





37 



Note. — Flux upper limits are given in units of 10 'cm ~s above 100 MeV and IGeV, for the 
power-law model eq.(T). 
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Energy (MeV) Energy (MeV) 

Fig. 1 . — Spectral fits to the counts (left panels) and the corresponding residuals (right panels) for the ROIs 
around two dwarf spheroidal galaxies, Willman 1 (top panels) and Draco (bottom panels). The lines in the 
spectral plots (left panels) are point sources (black), the Galactic diffuse component (blue) and the isotropic 
component (red). The black line overlaid to the data points is the best-fit total spectrum in the respective 
ROIs. The best-fit power-law models (with T = 2 here) for the dwarfs are below the lower bound of the 
ordinates. Willman 1 is the worst residual obtained in our sample, while Draco is illustrative of the fit quality 
for most ROIs. 
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3. Dark Matter constraints from dSph observations witli tlie Fermi-hAT detector 
3.1. Modeling of the dark matter density profiles 



While power-law spectra can be justified on astrophysical grounds, a proper search for a Dark Matter 
signal should take account of the specific spectrum resulting from WIMP annihilations. At a given photon 
energy E, the y-ray flux originating from WIMP particle annihilations with a mass niwiMP can be factorized 
into two contributions (IBaltz et al.ll2008b : the "astrophysical factor" J{tfr) related to the density distribu- 
tion in the emission region and the "particle physics factor" O^^ which depends on the candidate particle 
characteristics : 

(PwiMpiE, «A) - m X '^''''{E) , (2) 

where ifr is the a ngle between the d irection of observation and the dSph center (as given in Table [TJ. Follow- 
ing notations of 



BaltzetaL 



tooA J{il/) and O^^ are defined as 



1 < crv > 



2 An 



m 



WIMP f 



y dNf 
^ dE 



(3) 



and 



Jl.O.i 



dmp\im, 



(4) 



where < crv > is the relative velocity times the annihilation cross-section of the two dark matter particles, 
averaged over their velocity distribution, and the sum runs over all possible pair annihilation final states /, 
with dNf/dE and Bf the corresponding photon spectrum and branching ratio, respectively. The integral in 
Eq.([4ll is computed along the line of sight (l.o.s) in the direction if/, and the integrand p{l) is the assumed 
mass density of dark matter in the dSph. 



For each galaxy, we model the dark matter distribution with a Navarro-Frenk- White (NFW) (jNavarro et al 



19971 ) density p rofile within the tidal radius, as is reasonab le for cold dark matter sub-halos in Milky Way- 
type host halos Jpiemand et allboOTI : Ispringel et allboosl) : 







for r < rt 
for r > rt 



(5) 



where ps is the characteristic density, is the scale radius, and is the tidal radius. The line-of-sight 
integral in Eq.© may be computed once p.,, and are known. The tidal radius of the dwarf's dark 
matte r halo in the gravitational p otential of the Milky Way is self-consistently computed from the Jacobi 
limit (IBinnev & Tremainelll987r) for each set of and values assuming a mass profile for the Milky 
Way given by Catena & Ulhol ( 2009 ). The sharp cut-off in the density profile is rather extreme, but it is 
conservative in the sense that it truncates the probability distribution of expected y-ray fluxes at the high 
end. 



The halo parameter s (p,„ rQ, and the re sulting J factor from Eq.(|4l), are estimated following the 
methodology outlined in [Martinez et al.l (120091) . The observed line-of-sight (l.o.s.) stellar velocities are 
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well-d escribed by a Gaussian distribution (IMunoz et alJl2005L l2006l : IWalker et alj|200l l2009l : iGeha et al. 
2009bh and we include the dispersion arising from both the motion of the stars and the measurement errors 



as 



Strigari et all (120071) 



(=1 



: exp 



err . + cr^ .) 

los,i m,i' 



1 {Vj-uf 

Los,i m,i 



(6) 



where (v,) are the individual l.o.s. stellar velocity measurements and cr„, are the measurement errors on 
these velocities. The mean l.o.s. velocity of the dwarf galaxy is denoted by u. The full set of astrophysical 
parameters is £/ = ps, r^, u, and we discuss the two new parameters and p below. The theoretical 
l.o.s. dispersion, crios, is the projection of the 3D velocity dispersi on on the plane of the sky and this is 
determined using the Jeans equation (see lBinney & Tremainelll987l) once £/ is specified. T^t is the stellar 
mass to light ratio and it sets the mass of the baryons in these dwarf galaxies given the stellar luminosity. 
The velocity dispersion anisotropy is yS = 1 - cr^ /o"^, where cr, and cr,. are the tangential and radial velocity 
dispersion of the stars (measured with respect to the center of the dwarf galaxy). We assume that /? is 
constant for this analysis. The probability of the astrophysical parameters, £/ given a data set {v,} is obtained 
via Bayes' theorem: 'P(^|{v,)) oc P({vi }\£/)P(£/). The prio r probability, Pj^/), for the halo parameters. 



r .,py} is based on ACD M simulations (iDiemand et all 120071 : ISpringel et al.ll2008l) and described in detail 



Martinez et all (120091) . For we take the prior to be uniform between 0.5 and 5, and for /J the prior is 



m 



uniform between -1 and 1. 



The astrophysical factor J after marginalization over all the parameters in £/ for each dwarf galaxy 
within an angular region of diameter 1° is given in Table 01 The chosen 1° region for the calculation of J 
is a good match to the LAT PSF at energies of 1 - 2 GeV where most of the models under consideration 
are best constrained. At lower energies, the PSF is significantly larger, but beyond 1 ° the dwarf dark matter 
density has a negligible impact on the overall J computation, and at higher energies, the statistics with the 
current data are rather limited. Note that, due to their uncertain nature as true dark matter dominated dSphs 
or large uncertainties in their dark matter content, the Segue 2, Willman 1, and Bootes II dSphs have not 
been considered in this analysis. In addition, new stellar data on Segue 1 and Bootes II are being currently 
reduced and will be used in a forthcoming publication. We also exclude Ursa Major I, Hercules, and Leo 
IV, because their J values are smaller than those of the rest of the sample, yielding a final sample of 8 dSphs 
used for the dark matter constraints. 

In principle, annihilations in cold and dense substructure in the dwarf galaxy halo can increase J. 
However, previous studies have sh own that this boost d ue to annihilations in substructure is unlikely to be 



larger than a factor of few (see e.g. Martinez et al. 



2009h. Similarly, a boost in the annihilation cross-section 



in dwarfs due to a Sommerfeld enhancement (e.g. lArkani-Hamed et al.ll2009h . where the annihilation cross- 
section depends on the relative velocity of the particles, would increase the expected gamma-ray signal and 
improve our constraints. In order to be conservative, we have not included either of these effects. 
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3.2. Constraints on the annihilation cross-section 



Using Eq.©, the results given in Tabled and the DMFit package (IJeltema & Profumdl2008ar) as im- 
plemented in the Science Tools, 95% C.L. upper limits on photon fluxes (above 100 MeV) and on < crv > 
have been derived as a function of the WIMP mass, for each dSph and for specific annihilation channels. 
Our choices for the pair-annihilation final states are driven by theoretical as well as phenomenological con- 
siderations: a prototypical annihilation final state is into a quark-antiquark pair. The resulting y rays stem 
dominantly from the decay of neutral pions produced in the quark and antiquark hadronization chains, and 
do not crucially depend upon the specific quark flavor or mass; in fact, a very similar y ray spectrum is 
produced by the (typically loop-suppressed) gluon-gluon final state. Here, for illustration we consider the 
specific case of a hh final state: this choice is motivated by the case of supersymmetric dark matter (see 
Jungman et al.ll 19961. for a review). In supersymmetry with R-parity conservation, the prototypical WIMP 
dark matter candidate is the lightest neutralino, the mass eigenstate resulting from the superposition of the 
fermionic partners of the hypercharge and SU(2) neutral gauge bosons and of the two neutral Higgs bosons. 
At moderate to large values of tanyS, if the lightest neutralino is bino-like (i.e. if the U(l) hypercharge 
gauge eigenstate is almost aligned with the lightest neutralino mass eigenstate), dark matter dominantly 
pair-annihilates into hh. Another final state that is motivated by supersymmetric dark matter is into t^t~, 
that dominates in the case of a low-mass scalar superpartner of the r lepton, as is the case e.g. in the so- 
called co-annihilation region of minimal supergravity (mSUGRA). An intermediate case with a mixed hh 
and t^t" final state is also ubiquitous in supersymmetry, since those are the two heaviest down-type matter 
fermions in the Standard Model. The additional color factor and a larger value for the mass typically favor 
a relatively larger hh branching fraction. While the choice of the final states we consider is motivated here 
by supersymmetry, the results we find apply to generic WIMP models and not only to neutralinos. 

Fig. [2] shows the derived upper limits on the photon fluxes for all selected dwarfs and for various 
annihilation final states (respectively, 100% hh in the upper left panel, 100% t^t" in the upper right panel, 
and a mixed 80% hh -\- 20%t'^t^ final state in the lower left panel). The lower right plot of Fig.[2]illustrates, 
for the case of Ursa Minor, how upper limits on the y-ray fluxes change as a function of mass depending on 
the selected final state. Final states producing a hard y-ray spectrum, such as fi^fi^ and t'^'t^ result in the best 
upper limits, since they predict abundant photons fluxes at larger energies, where the diffuse background is 
lower. With increasing mass, the spectrum from WIMP annihilation is rigidly shifted to higher energies, 
and the advantage of having a harder spectrum is less critical: for niDu ~ 1 TeV the upper limits we obtain 
for all the final states we consider vary within a factor 3, versus more than one order of magnitude at lower 
masses. 

The results presented in fig. |2] bracket realistic cases of theoretically motivated particle dark matter 
models, but since they only refer to the WIMP annihilation final state they do not depend on the particular 
assumed particle theory model. We consider next a few motivated specific WIMP dark matter scenarios, 
and study hov y the corresponding pararneter space is constrai ned. We consider first the well-known case 
of mSUGRA dchamseddine et all (ll982l) : lBarbieri et all Jl982l) : see alsolNiUsa (11984) for an early review), 
where the supersymmetry breaking parameters are typically specified at a high-energy scale (typically taken 
to be the grand unification scale Mqut - 2 x lO'^ GeV) and assumed to be universal at that scale. Those 




Fig. 2. — Derived upper limits on fluxes for all selected dwarfs and for various branching ratios: 100% bb 
(upper left), 100% t^t" (upper right) and mixed 80% bb + 20%t^t^ (lower left) final state. Lower right 
plot gives an illustration of how the upper limits on the fluxes can change depending on the selected final 
state (here for the Ursa Minor dSph). 
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parameters are the universal high-scale supersymmetry-breaking sfermion mass mo, gaugino mass M1/2 
and trilinear scalar coupling Aq; an additional (low-scale) parameter is the ratio of the vacuum expectation 
values of the two Higgs doublets, tan/3, and the sign of the higgsino mass parameter //. Here, we adopt the 
following ranges for the various parameters (linear scan) : 80 < niQ < 100000, 80 < M1/2 < 4500, Aq = 0, 
1.5 < tan/3 < 60 and with sign(//) being not constrained. 

A less constrained alternat ive is to consider a "phenomenological" Minimal Supersymetric Standard 
Model (MSSM) setup (see e.g. IChung et al.ll2005l for a review of the most general MSSM Lagrangian), 
where all soft supersymmetry breaking parameters (i.e. the positive mass-dimension coefficients of la- 
grangian terms that explicitly break supersymmet ry) are defined at the electr o-weak (low-energy) scale, pos- 
sibly with a few simplifying assumptions (see e.g. Profumo & Yagunall2004 , for an early atter npt at a scan of 



the MSSM parameter space). Here, we consider the reduced set of parameters considered in Gondolo et al. 



and perform a logarithmic scan over the following parameters : < 10000, IM2I < 10000, 
100 < rriA < 1000, 1.001 < tanyS < 60, 100 < nig < 20000, while the scan is linear in -5 < A,/mg < 5 and 
-5 < At/m^ < 5. 

In addition to the two scenarios considered above, we also entertain two additional specific WIMP dark 
matter m odels. The first on e is the lightest Kaluza -Klein particle of Universal Extra Dimension s (UED) 
(see e.g. ICheng et al.1 (|2002h : IServant & Taiti (|2003h . for a review see iHooper & Profumol (120071) ). where, 
in the minimal setup, the dark matter candidate corresponds to the first Kaluza-Klein excitation of the U(l) 
hypercharge gauge boson, also known as B^'^ In this case, there is an almost fixed relationship between the 
dark matter mass and its pair annihilation cross section, an d a thermal relic abun dance in accord with the 
dark matter density is obtained for masses around 700 GeV (|Servant & Taitil2003h . 



The second model, which was recently considered in Kane et al.l (l2009h as a nat ural and well motiva ted 
scenario in connection with the anomalo us positron fraction repo rted by PAMELA (lAdriani et al.ll2009h . is 
that of wino-like neutralino dark matter (|Moroi & Randallll2000l) . Wino-like neutralinos (which for brevity 
we will refer to as "winos"), i.e. neutralino mass eigenstates dominated by the component corresponding 
to the supersymmetric fermionic partners of the SU(2) gauge bosons of the Standard Model, pair annihilate 
very efficiently into pairs of W~ (if their mass is larger than the W mass) and the pair annihilation cross 
section is fixed by gauge invariance once the wino mass is given. Winos arise in various supersymmetry 
breaking scenarios and in several string motivated setups, where e.g. the anomaly mediation contribution 
to the gaugino masses dominates over other contributions, setting M2 < Mi . Typical such scen a rios o f 
anomaly med i ated s upersymmetry breaking (AMSB) were considered e.g. in lRandall & SundrumI (Il999h : 
Giudice et all (Il998h . Although winos produce a thermal relic density matching the universal dark matter 
density for masses around 2 TeV, several non-thermal production mechanisms have been envisioned that 
could explain a wino dark matter scenario with lighter dark matter candidates than a TeV. 

Fig. |3] compares the resulting LAT sensitivity in the (mwimp,< crv >) plane with predictions from 
mSUGRA, MSSM, Kaluza-Klein dark matter in UED and wino-like dark matter in AMSB. All mSUGRA 
and MSSM plotted models are consistent with accelerator constraints. Red points are compatible with the 
3cr WMAP constraint on the universal matter density under thermal production while blue points would have 



-17- 



a lower thermal relic density. For the blue points, we assume that the production mechanism is non-thermal 
in order to produce the observed universal matter density, and we therefore do not rescale the neutralino 
density by the factor (Hthermai/f^DM)^. which would result assuming exclusive thermal production. This is 
very natural in the context of several string-theory motivated frameworks, where moduli generically decay 
into both Standard Mo del particles and their sup ersymmetric partners, which in turn eventually decay into 
the lightest neutralinos Moroi & Randalll (|200(J) . Topological objects such as Q-b a lls ca n als o decay and 
produc e neutralinos out of equilibrium, as envisioned e.g. by iFujii & Hamaguchil (|2002h and iFujii & Ibe 
(|2004) . Another possible scenario is one where the expansion history of the Universe is more rapid than in a 
radiati on dominated setup, for instance bec ause of a dynamical "Quintessence" field in a kinetic-dominated 
phase JsalatilbooilProfumo & UUiolboOBh . 

Fig. |3] clearly shows that, after less than a year of Fermi data survey, the upper limits on the y-ray flux 
from dSphs are already starting to be competitive for MSSM models, provided that these models correspond 
to low thermal relic density. Draco and Ursa Minor dSphs set the best limits so far. Pending more data, they 
may also start to constrain mSUGRA models with low thermal relic density as well. Furthermore, these 
flux upper limits already disfavor AMSB models wit h masses <300 GeV. Interestingly, our results strongly 
constrain the models considered in lKane et al.l (l2009f) . invoking a 200 GeV mass wino. 



3.3. Comparison to dark matter models proposed to fit the PAMELA and Fermi e+e data 



The recent detection by the PA MELA experiment of a positron fraction that increases with energy 
above 10 GeV (|Adriani et al.l 120091) and the possibility that dark matter annihilation in the Galaxy could 
produce this "positron excess" (among more mundane explanations such as pulsars) have spurred great 
interest in the particle physics community. The pair annihilation of galactic WIMP dark matter can, in 
principle, produce an anomalous excess in the positron fraction at energies between a few GeV and -100 
GeV The spectrurn of high-energy + e~ , although compatible with a purely canonical cosmic-ray origin 
(lAbdo et al.ll2009l : iGrasso et al.ll2009r) . can also accommodate an additional component due to galactic dark 
matter annihilation. 

A dark matter annihilation interpretation of the positron excess implies preferentially a leptonic final 
state, to avoid the over-production of antiprotons. This is very hard t o achieve in the minimal supersymmetric 



Kane et al. 



20091. for the case of AMSB, which we will 



extension of the Standard Model (see however 
not consider further here). In addition, the spectral shape of high-energy e'^e^ points towards rather large 
masses, and the level of the needed local positron flux indicates either a very large pair-annihilation rate, or a 
strong enhance ment in the local d ark rn atter density. Using a canonical primary electron injection spectrum, 
the analysis of iBergstrom et al.l (120091) further indicates that a preferred annihilation final state is fi^fi', 
or the somewhat softer (in the produced e'^e^ spectrum) but essentially very similar four body fi^ fi^ fx^ fi~ 
final state. Theoretical arguments that could explain this peculiar annihilation fin al state, possibly invo lving 
mechanisms to enhance the low-velocity annihilation rate, have been proposed (Arkani-Hamed et al.ll2009, : 
Nomura & Thaleiil2009r) . With standard assumptions on the dark matter density profile, and assuming a 
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Fig. 3. — mSUGRA (upper left), MSSM (upper right), Kaluza-Klein UED (lower left) and Anomaly me- 
diated (lower right) models in the (mwimp,< crv >) plane. All mSUGRA and MSSM plotted models are 
consistent with all accelerator constraints and red points have a neutralino thermal relic abundance corre- 
sponding to the inferred cosmological dark matter density (blue points have a lower thermal relic density, 
and we assume that neutralinos still comprise all of the dark matter in virtue of additional non-thermal pro- 
duction processes). The lines indicate the Fermi 95% upper limits obtained from likelihood analysis on the 
selected dwarfs given in Table ID 



Table 4. Properties of the dark matter halos of dwarf spheroidals used in this study. 



Name 




[<«2> - <i?)2, <p2) - <p>2, (RP) - (R){P) 
R = log,o(r,//:pc), P = log,o(p.,/Mo kpc' 


1 jNFW 
3) (lOl'Cfld) 

Clip 


Ursa Major II 
Coma Berenices 
Bootes I 
Usra Minor 
Sculptor 
Draco 
Sextans 
Fornax 


[-0.78, 8.54] 
[-0.79, 8.41] 
[-0.57, 8.31] 
[-0.19, 7.99] 
[-0.021, 7.57] 
[0.32, 7.41] 
[-0.43, 7.93] 
[-0.24, 7.82] 


[0.0417, 0.0986, -0.0554] 
[0.0603,0.132,-0.0820] 
[0.0684, 0.165, -0.0931] 
[0.0430, 0.116, -0.0697] 
[0.0357, 0.0798, -0.0528] 
[0.0236, 0.0364, -0.0286] 
[0.0302, 0.109, -0.0570] 
[0.0474, 0.140, -0.0798] 
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Note. — These parameters are obtained from measured stellar (line of sight) velocities, and Tj are 
the density and scale radius for the dark matter halo distribution. The first column, [logmlpj), logjoCr,)], 
is the average in the joint logio(r,) — logmtp,) parameter space, whose posterior is well described by a 
Gaussian distribution centered on the average value given. The second column gives the diagonal and 
off diagonal components of the covariance matrix that may be used to approximate the joint probability 
distiibution of p, and as a Gaussian in log[Q(r,) and logio(p,). The last column provides J^^^ (see 
Eq.|4), which is proportional to the pair annihilation flux coming from a cone of solid angle 2.4 10""* sr 
centered on the dwarf. The errors on J^^^ are obtained from the full MCMC nrobabilitv distribution 




Fig. 4. — Constraints on the annihilation cross-section for a /j.'^'/j." final state based on the 95% confidence 
limits on the y-ray flux compared to dark matter a nnihilation model s which fit well either the PAMELA 
(lAdriani et al.ll2009 ) or Fermi + e~ measurements (lAbdo et al. The left panel shows the constraints 

considering y-ray emission from final state radiation only. The right panel shows the constraints for the Ursa 
Minor dwarf including both y-ray emission from IC scattering and final state radiation. Here we consider 
two different diffusion coefficients, and show the effect of the uncertainties in the Ursa Minor density profile. 
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final state, the regions in the pair-annihilation cross section versus mass plane preferred by the Fermi- 
LAT e^e' data are shown in orange in Fig.|4l while those favored by the PA MELA positron fraction dat a are 
highlighted in light blue (for details on the computation of these regions see lProfumo & Jeltemal (l2009h ). 



In a pair-annihilation event producing a ji^if pair, y rays result from both the internal bremsstrahlung 
off of the muons (final state radiation), with the well-known hard power-law spectrum AnylAEy ~ E~\ and 
from the inverse Compton (IC) up-scattering of cosmic microwave background light by the e^e' resulting 
from muon decay. The dark matter interpretation of the "cosmic-ray lepton anomalies" implies significant 
y-ray emission from a variety of sources; pr edictions and constraints o n these models have been discussed 
extensively in the recent literature (see e.g. IProfumo & Jeltemal l2009l. for a discussion of the constraints 
from the expected IC emission from annihilation at all redshifts and in all halos). 

The left panel of Fig. |4] illustrates the constraints we derive from 1 1 months of Fermi data on local 
dSph on a generic WIMP dark matter pair-annihilating into a jj^fi' final state (we do not specify here any 
particular particle physics scenario, although as stated above several examples have been considered in the 
literature), considering the final state radiation emission only. Here the spectral modeling was done using 
the DMFit package as in 93.21 The hierarchy among the constraints derived from the various dSph is very 
similar to what we find for the softer pair-annihi lation final states consider ed in Fig. |3] Neglecting any 
low-velocity enhancement of the annihilation rate (|Arkani-Hamed et al which would boost the y-ray 

signal from dSph and hence the constraints we show, the final state radiation alone does not yet exclude 
portions of the parameter space favored by the dark matter annihilation interpretation of the cosmic ray 
lepton data. 

The calculation of the y-ray yield from IC is complicated, in systems as small as dSphs, by the fact that 
e'*'e~ are not confined (i.e. their diffusion lengths are typically larger than the physical size of the system). 
In fact, the typical energy-loss length-scale for TeV e'*'e' loosing energy dominantly via IC off of photons 
in the cosmic microwave background is of the order of hundreds of kpc, much l arger than the size o f dSph . 
Assur nptions on cosmic ray diffus ion are therefore critical, as discussed e.g. in IColafrancesco et all (|2007h 
and in ljeltema & Profumd (|2008br) . In the absence of any direct cosmic-ray data for external galaxies such as 
dSph (the only piece of information being that dSph are gas-poor environments with typically low magnetic 
fields), we consider the usual diffusion-loss equation and solve it in a spherically symmetric diffusive region 
with free-escape boundary conditions. We employ a diffusion coefficie nt at the level of what is usually 
inferred for cosmic rays in our own Milky Way (e.g. IStrong et al.ll2007h (for a th oroug h discussion of the 
diffusio n model we adopt here, we refer the reader to Colafrancesco et al. ( 200?!) and jjeltema & Profumo 
(l2008bl) ). Specifically, we consider two values for the diffusion coefficient. Do = 10^^ cm^/s and Dq = 
1029 cm /s bracketing the values typically inferred for the Galaxy: the larger the diffusion coefficient, the 
larger the cosmic ray mean free path, and the larger the leakage of cosmic ray e'^e~ out of the dwarf, leading 
to a suppression of the IC signal. We also assume a power-law dependence of the diffusion coefficient on 

I In the energy loss term, IC emission off of Cosmic Microwave 



energy given by D{E) = Dq (^-^ 



, 1 GeV / 

Background photons by far dominate over both synchrotron and starhght IC losses. 



We show our results, including both final state radiation (FSR) and IC emission offofCMB photons, 
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in Fig.m right panel, for the case of the Ursa Minor dSph. Here the spectrum was modeled self-consistently 
with custom spectra which include the expected y-ray emission from both IC and FSR for a given assumed 
particle mass and diffusion coefficient for a grid of particle mass values ranging from 100 MeV to 10 TeV. 
The Ursa Minor dSph was chosen as an example of one of the best cases for this particular study. In the 
smaller ultra-faint dwarfs, diffusion is expected to have a much larger effect due to the much smaller diffusive 
region (modeled based on the stellar extent), and IC emission for the diffusion coefficients assumed here will 
not add significant flux above what is expected from FSR alone. 

The lower and upper lines in Fig. |4] indicate the range of uncertainty in the determination of the dark 
matter density profile of the Ursa Minor dSph. With the smaller diffusion coefficient choice, and for large 
enough masses (producing higher energy e'^e~ and subsequent IC photons), the IC emission dominates, and 
it exceeds our y ray upper limits for models that fit the PAMELA data and have masses larger than 1 TeV. 
Excessive emission is predicted also in the more conservative case with Dq = 10^^ cm^/s for some models 
with masses in the 2-5 TeV range. 



4. Conclusions and final remarks 



We have reported the observations of y-ray emission from 14 known dwarf spheroidal galaxies by 
/^erm/-LAT. No excesses have been observed in LAT data and upper limits have been derived on the y-ray 
flux from dSphs. 

Using the dark matter halo modeling for the 8 best candidate dwarf spheroidal galaxies derived from the 
latest stellar data (tab.©, we have shown that if dark matter is assumed to consist entirely of neutralinos, the 
upper limits obtained from one year of LAT data begin to constrain mSUGRA and MSSM models with low 
thermal relic densities and AMSB models with wino-like neutralinos with masses below 300 GeV (fig.O. It 
is worth noting that fou r dSphs have also been observed by Cherenkov telescop es : Sagittari us by H.E.S.S . 
(|Aharonian et allboOSt) . Draco and Ursa Minor by Whipple dwood et allboosb and Ve ritas Jwagneil 12009 . 



which also includes Willman 1), and finally Draco (|Albert et al 



2008h and Willman 1 ( Aliuetal 



I200^ bv 



MAGIC. The observation time varies between these studies, but in general the limits on the annihilation 
cross-section reported vary between a few times 10"^^ to a few times lO"'^^ cm~^ s~^ for a 1 TeV mass neu- 
tralino and an assumed NEW dwarf density profile. lACT observations are most sensitive to typically higher 
mass dark matter particles (greater than ~ 200 GeV) compared to the LAT, making them complimentary to 
Fermi searches. 

The Fermi limits also constrain WIMP models proposed to explain the Fermi and PAMELA data, 
particularly for high particle masses (> 1 TeV, fig. HJl. For these models, strong constraints come from the 
inclusion of the expected IC y-ray emission, though the flux of this component depends on the assumed 
diffusion model of e'^e~ in dSphs. 

It is worth emphasizing that the results presented in this paper have all been obtained for a standard 
NEW halo shape and without assuming any boost factor effect due to substructures in the dwarfs or a 
Sommerfeld enhancement to the annihilation cross-section. 
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